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We study the influence of a chaotic environment in the evolution of an open quantum system. We show that 
there is an inverse relation between chaos and non-Markovianity. In particular, we remark on the deep relation 
of the short time non-Markovian behavior with the revivals of the average fidelity amplitude - a fundamental 
quantity used to measure sensitivity to perturbations, and to identify quantum chaos. The long time behavior is 
established as a finite size effect which vanishes for large enough environments. 
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I. INTRODUCTION 

The advent of quantum information and quantum technol- 
ogy has brought us to a deeper understanding of the basics of 
quantum mechanics. As a result, various fundamental aspects, 
such as equilibration [1], simulability |2|], and even founda- 
tions (Hi, have been revised. This framework has provided 
both valuable insight into the foundations and new techno- 
logical achievements. But theory and its related experiments 
have advanced asymmetrically mainly due to the impossibil- 
ity to isolate completely the experimental setup, leaving the 
system 'open' and exposed to decoherence fl]. 

Many quantum open systems problems can be solved after 
assumptions are made. The widely used Born-Markov ap- 
proximation, successfully applied to describe many physical 
situations J5[], is associated with both a memory-less envi- 
ronment and a weak coupling between system and environ- 
ment. However, recently, interest in quantum open systems 
where this assumption no longer applies - usually called non- 
Markovian (NM) evolution - has flourished ]6|48|] . 

A natural question to ask is to what extent the dynamical 
properties of the environment can extend the validity of the 
Born-Markov approximation, even beyond the weak coupling 
regime. In other words, how well does a chaotic environment 
reproduce Markovian evolution. In this paper, we address this 
question analytically and numerically by means of a probe 
qubit coupled to a generic environment where different de- 
grees of chaos can be tested. Exploring the time evolution 
of non-Markovianity measures, we show that the stronger the 
chaos in the environment, the more Markovian the evolution 
is. We build upon previous knowledge [|9l|l0t] of measurable 
quantities such as the fidelity amplitude which can be directly 
related to recently proposed measures of NM behavior. More- 
over, we establish that the short time behavior of the fidelity 
amplitude determines the characteristics of the qualitative NM 
behavior. The lingering, long time contributions, are finite 
size effects which contribute only as a linear term with a slope 
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that goes to zero as the environment size goes to infinity. The 
remaining non-Markovianity is thus size independent and due 
to short time revivals in the fidelity amplitude decay. 

The paper is organized as follows. In Sect. [II] we describe 
the system we use for our analysis. It consists of a qubit in 
the presence of an environment whose evolution is subject to 
the state of the qubit. In Sect. [Ill] we briefly describe the Non- 
Markovianity measure based on distinguishability and its re- 
lation to the fidelity decay of the environment. The numerical 
results and analysis are done in Sect. [IV]and we present some 
concluding remarks in Sect. |V] 



H. SYSTEM 

We consider a system-environment situation, whose Hilbert 
space is ,W — (£> ^ n v, where sys and env denote system 
and environment. Any Hamiltonian can then be split as 

" = ^sys ~\~ Hem ~\~ £^sys,env ; ( 1 ) 

where £ is real and controls the strength of the interaction be- 
tween system and environment. We further restrict to the case 
in which Jf sys represents a qubit, i.e. when dimJf sys = 2, and 
we call N = dimJff env . The choice of the terms in Eq. (fl} 
results in different physical situations. However, instead of 
specifying the particular form of Eq. (Q3, we shall impose two 
general conditions. The first one involves the nature of the 
interaction. We shall assume that it is factorizable, i.e. that 
Vsys,env = Vsys <8> V e nv- Such structure appears in a wide va- 
riety of situations, including Ising interaction and atom-field 
interaction under various approximations ill ill . The second 
assumption is to consider that the evolution of the central sys- 
tem, Jtlys, either occurs at much smaller time scales than that 
in which decoherence occurs or that it is a multiple of V sys . 
In the former, one can safely ignore the contribution to the 
dynamics, and in the latter case, H sys can be included in the 
interaction term and can keep the factorizable structure un- 
affected. This occurs, e.g. in the case of a strong magnetic 
field applied to a set of interacting spins and thus is of particu- 
lar importance in, among others, nuclear magnetic resonance 
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(NMR). In this situation, one can write 

H=|0>(0|®flb+|l)<l|®£ri 



(2) 



with Hq and Hi acting only on the environment and both 
0}(0|, being projectors onto some orthonormal basis of 
the qubit. This Hamiltonian has already been introduced in 
lfl2tl and can be interpreted as having an environment whose 
evolution is conditioned by the state of the qubit. The initial 
state of the system shall be p sys ,env (0) = p S ys (0) ® p cm - Notice 
that the only condition imposed on the initial state is that it be 
a product state. The evolution of the qubit is 



with 



Psys(0 = Trenv [U(t) Psys (0) ® Pem UHt)] (3) 



U(t) = \Q)(0\U (t) + \l)(l\U 1 (t) (4) 



and {//(f) =exp(—itHj/h). It is convenient to rephrase Eq. (0 
in terms of a quantum channel 



p sys (0=A(*)[P S ys(0)] 



(5) 



Notice that although A(f ) results from tracing out the environ- 
ment it still depends on p env . 

The matrix elements of the channel induced by Eq. (01, in 
the Pauli basis are 



A S = i Tr [0}U{f)a k ® Penvf/ + (f )] 



V.* 

If we take Co = I and 01,2,3 
simple form 



(6) 



Ox,y,z> the channel takes the 



A = 



no o\ 

9t[f(f)] 3 [/(f)] 

3 [/(f)] 9t[/(r)] 

\0 1/ 



(7) 



with /(f) = Ti[p em Ui(tyUo(t)] being the expectation value 
of the echo operator. If p env is a pure state [10] then |/(f)| 2 
is the Loschmidt echo (LE) lfl3ll - also called fidelity - orig- 
inally proposed to measure sensitivity to perturbations in the 
Hamiltonian as a signature of quantum chaos iffl . The LE 
decays as a function of time and the - more or less - univer- 
sal decay regimes have been extensively studied (see reviews 
lfl5l,[l6t0 . The environment could in fact be in a pure state - 
e.g. in a thermal ground state at zero temperature. However 
it is probably easier to imagine the environment being in a 
mixed state - e.g. at thermal equilibrium at a given tempera- 
ture. We choose then the environment to be in the maximally 
mixed state, i.e. proportional to the identity. In that case we 
obtain the real and imaginary part of the average fidelity am- 
plitude (AFA), 



(/(f)) = -Tr[f/ 1 (f)tf/ (f) 



(8) 



i.e. the average value of the echo operator with respect to an 
orthonormal basis. We remark that the choice of basis (or any 
complete set) is arbitrary. This fact contrasts the case of the 
LE where the kind of states in the set is crucial 12111 . 



III. FIDELITY AMPLITUDE MEASURES 
NON-MARKOVIAN BEHAVIOR 

During a classical Markovian process the distance between 
two initial distributions decreases monotonically. Devia- 
tions from this behavior are a landmark of non-Markovianity. 
Breuer et al. [7] used this property to define a measure of NM 
behavior in a quantum setting. The distance can be chosen as 
to link non-Markovianity with distinguishability of states and 
thus information flow between the system and its surround- 
ings. Such a measure is defined as 



max / dta{t,p\ 2(0)), 

Pl,2(0)i(T>0 



(9) 



where ff(f,p 1-2 (0)) = dD(p l (t),p 2 (t))/dt is the rate of 
change of the trace distance 

D(p 1 (f),p 2 (f)) = itr|p 1 (f)-p 2 (f)| (10) 

between initial states p 1,2(0). In ||8[] two other measures 
were proposed, based on deviation of semi-group properties 
of quantum flows. Both study the physicality of the induced 
instant map at intermediate times, one via the Jamiolkowski 
isomorphism and the other via the entanglement (as measured 
with the concurrence) with an ancilla qubit. It is straightfor- 
ward to show that for channels like Eq. (0 the measure in- 
duced by the entanglement is proportional to „# . 

In our case, the states that maximize ./# are p± = (I ± 
(aa x + ba y ))/2, with |a| 2 + \b\ 2 = 1, leading to 



= 2 



l/l>o 



d\(f(t))\ 
dt 



(ID 



In other words, Eq. (fill means that NM behavior is directly 
related to the positive derivative of the AFA as a function of 
time lfl7tl . 

For fully chaotic systems both the AFA and the LE saturate 
around a value that depends on h. After saturation sets in, the 
state is approximately random and the value of both fidelity 
and fidelity amplitude fluctuate. As a consequence, we expect 
the NM measure to grow indefinitely. Thus in our calculations 
of ^# we modify the original definition in Q and calculate 
the NM measure up to a certain time. Regardless, the measure 
at time f still holds its meaning, i.e. the larger ^#(f ) means 
the distance between the two states has ceased to decrease (or 
increased) more in that period of time, which implies a more 
NM behavior. 



IV. NON-MARKOVIANITY AND CHAOS: RESULTS 

Now we consider the long standing question of the relation 
between chaos and Markovianity. To do so we model the en- 
vironment using simple but fully featured systems: quantum 
maps on the torus. The quantization of the torus implies that 
both position and momentum are discretized and the effective 
Planck constant is the inverse of the Hilbert space dimension 
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FIG. 1. (color online) JV it) for the PCM (top) and the HM (bottom). 
Top: (squares) A = 0.96, (circles)A = 1.76, (triangles) A = 5.99. 
5K/Ti = 3.635. (diamonds) A = 0.96, 5K jh = 0.64 (weak coupling). 
Bottom: (squares) k = 0.001 (regular), (circles) k = 0.25 (mixed), 
(triangle) k = 1 (fully chaotic). Top inset: (diamonds) weak coupling 
case. Botom inset: (triangles) k= 1, note the linear dependence. The 
slopes of the straight lines are a = 1/2N (top panel, dashed) and 
a = 2/3N (bottom inset, dot-dash). N = 4096. 



N. In this setting, a quantum map is simply a unitary U acting 
on an N dimensional Hilbert space. 

We consider two different maps. First, the quantum per- 
turbed cat map (PCM) 



Ur K = 



—iKNap ixNaq-ixNK(2sm(2nq) — s,m(A7iq)) 



(12) 



where q and p are the generators of periodic position and mo- 
mentum translations on the torus with discrete eigenvalues 0, 
\/N, ...,(N -\)/N. The subindex K denotes the depth of the 
kicking potential. For a = 1 and K = it is the quantum ver- 
sion of Arnold's cat map, a uniformly hyperbolic and mixing 
map of the torus onto itself, which is a paradigmatic exam- 
ple of chaos in two dimensions. The positive Lyapunov expo- 
nent A, which determines the rate of exponential divergence of 
classical trajectories, is uniform over the whole phase space. 
We explore different degrees of chaos by changing a since, for 
small K, A « ln((2 + a 2 + y/a 2 (4 + a 2 ))/2). Here 17(f) = U' 
where now t is an integer, and Uq = U c .k and U \ = U c k+Sk- 
The other map we consider is 



^HJfcJf 



JNkcos(2nq) iNk 1 cos(27rp) 



(13) 



which corresponds to the Harper map (HM) [18]. It is an ap- 
proximation of the motion of an electron in a crystal under the 
action of an external field 11911 . For A: < 0.11, the dynamics 
described by the associated classical map is regular, while for 
k > 0.63 there are no remaining visible regular islands. We 
set Uq = U H ,k,k and U\ = U HMSk ,k- 

We now take the result of Eq. (fTTT i and compute numeri- 
cally ^#(/) for the two maps. Notice that the structure of the 
maps is U\ = UqP(e) (with e = 8K or 8k ). The e — > limit 
implies P(e)[and/(f)] — > 1 (i.e. no decoherence). The cou- 
pling strength is given by 8K and 8k. For weak couplings, 
(in the chaotic case) the AFA decays exponentially, and the 



O o nu DDDn 




ooooooooooooooo 



o o 



□ □ 




10 15 20 25 30 35 40 
t 



FIG. 2. (Color online) |(/(f))|, same cases as Fig. 1 for the PCM 
(top) and the HM (bottom). Inset: A = 0.96 for the PCM, for N = 
512,4096, 16384. 



rate depends quadratically on the coupling parameter - the 
Fermi golden rule regime (FGR). Here the evolution is ex- 
pected to be Markovian. Throughout this paper we consider 
mainly coupling strengths beyond the FGR. 

The AFA is evaluated directly by averaging the echo oper- 
ator over a complete set of states. For the PCM we change 
a in Eq. ( fT2l so we can assess the change in ^# for differ- 
ent levels of chaoticity. For the HM varying k in Eq. ( fl"3l ). 
we go from integrable to completely chaotic. In Fig. Q] (top), 
we show the NM measure ^#(f) for three different examples 
of the PCM with varying degrees of chaoticity. On the bot- 
tom of Fig.Q] we show the same for the HM, where we show 
the AFA for three different k (regular, mixed, and chaotic). 
For the PCM after a small number of steps, there appear three 
distinct jumps. As expected, the larger A is, the smaller is 
the jump, which confirms the intuitive relation that the more 
chaotic the environment is, the more Markovian is the evo- 
lution. After this short time behavior, the three cases ex- 
hibit linear growth of j$t{f). The explanation is simple. For 
fully chaotic systems, at a time of the order of Ehrenfest 
time [?e = ln(l//i)/A], the AFA saturates - but oscillates - 
around l/N This saturation corresponds to the overlap be- 
tween two completely random states, and is approximately 
constant. This implies that will grow linearly and that the 
slope a will be proportional to 1 /N. The proportionality con- 
stant depends on the map (for the PCM, regardless of the value 
of A, we found the slope to be a w 1 /2N). For completeness, 
we include the case A = 0.96 in the weak coupling regime. 
As expected, ) = visibly up to t sa 50. After that we 
have linear growth. On the bottom of Fig. Q] we show j$(f) 
for the HM for three qualitatively different cases. In the case 
in which the classical dynamics is regular (k = 0.001), we ob- 
serve that up to the times shown j${t) increases non-linearly. 
There is a saturation of the AFA, but not at 1 /N, which even- 
tually leads to a linear growth of jfc '. This saturation for small 
k takes place at much larger times. When the dynamics is al- 
most fully chaotic (k = 1), there is a very small jump after 
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which there remains only the linear growth due to fluctuations 
around the saturation value. The slope of this linear growth 
is a « 2/ 3N. In the parameter region where the KAM tori of 
the HM begin to break, there is a combination between regu- 
lar and chaotic dynamics (initial states can have components 
inside regular islands and components inside the chaotic sea) 
and the behavior is less intuitive. In fact what is observed - 
in Fig. Q]for k = 0.25 -is that the environment modeled by a 
HM in the transition from regular to chaotic can be strongly 
non-Markovian [see also Fig.[3). 

In both situations, the long time behavior for the NM mea- 
sure is linear. This would imply ^# — > °°. However, this as- 
sertion presents no problems in our analysis. The slope of 
the long time linear regime goes to zero as N grows. Intu- 
ition suggests 'large' environment as a necessary condition of 
Markovianity. However, in the N — > °° limit, there will always 
remain the short time value attained by (see Fig.Q~|i, which 
is independent of N (Fig. |2). 

To shed more light on the results displayed in Fig. Q]we fo- 
cus on the evolution of the AFA as a function of time for the 
cases considered above. In Fig. [2] we show examples of the 
decay of the square of the AFA. In the top panel, we show 
results for the PCM for three different values of X (i.e. differ- 
ent a in Eq. (fT2b). In the bottom panel, we show the same for 
the HM, with three values of k (regular, mixed, and chaotic). 
In contrast with the LE, the AFA is independent of the type 
if initial states and decays exponentially with two distinct de- 
cay rates. The short time decay rate Y can be related to un- 
corrected - random - dynamics ||9|. The value of Y can be 
computed using semiclassical methods. This decay rate can 
diverge, meaning that the short time decay can be extremely 
fast. These divergences - which depend on the type of per- 
turbation, and are more evident the larger X is - could be re- 
lated to the phenomenon known as survival collapse l20ll after 
which the largest revivals appear. For the numerical results, 
on the top we chose a value of 8k for the PCM which cor- 
responds to a large Y (near the div erging values), where the 
largest revivals have been observed ll9ll2lll. We remark, more- 
over, that the short time decay of the AFA is independent of 
N and therefore so is the revival. In the inset of Fig. |2]we see 
that the AFA (for the PCM, with X = 0.96) is almost equal 
for three different values of N up to t « 10. This is important 
because the short time revivals will provide the main contri- 
bution to ./#. While this contribution remains constant with 
N the long time contribution goes to zero as \/N. The curve 
with diamonds supports the results shown in Fig. Q] [top] for 
the weak coupling regime. 

The possibility to assess the behavior as an environment 
model, by changing one parameter, from regular to chaotic, is 
indeed tempting. In Fig. [3] we computed ^# at a fixed time 
for the HM for different values of k. We chose t — 20, around 
the time in which the fastest decaying case starts to saturate 
(see Fig [2] bottom). We see that for small k (regular dynam- 
ics), ^ takes a constant value (which, apart from the fixed 
time, depends on N and 8K) and there is a transition where 
^ depends on k just where the KAM tori begin to break. 
When the dynamics is fully chaotic, the value of ^# (at t = 20) 
again takes a constant value. Fig. [3] is a clear example of the 



k = 0.001 k = 0.25 k = 1 



0.6 
M 0.4 



0.1 











• \ 

I 


1 

1 — , , — 1 — 1 — , , — 1 — ' — . . — 1 — ' — ' "0 "0 — e— e ©-J 



1 — ' — ' — ■ — ' "c "o — e— ©- 

10- 4 to -3 io~ 2 10- 1 io n 10 1 



k 

FIG. 3. Jt{t = 20) as a function of k, for t = 20, for the HM with 
N = 4096 and Sk — 0.00113. Top: phase space diagram for three 
examples of k: regular, k = 0.001; mixed k = 0.25; chaotic k= 1. The 
corresponding points in the bottom panel are drawn in solid black. 



expected behavior: regular environments are expected to be 
more NM while the NM behavior that appears to linger in the 
chaotic regime is due to the same oscillatory behavior around 
the saturation value mentioned for the case of the PCM. In 
the transition region, 0.11 < k < 0.63, there is coexistence 
between tori and chaotic regions. In the first place, the exis- 
tence of regular islands implies that even though there will be 
leaking - by tunneling - to the chaotic regions, the saturation 
will take much longer. In addition, the area occupied by the 
chaotic region is smaller than the torus, therefore, the satu- 
ration value is larger than 1 /N. We have checked for other 
times (up to t ~ 1000) and also other methods (not shown) 
- e.g. taking as NM value the y-intercept of an asymptotic 
linear fit - and the qualitative behavior is the same. Further 
studies are needed in order to fully grasp the behavior in the 
intermediate region. 



V. CONCLUSIONS 

We addressed the issue of how well a chaotic environment 
can model Markovian evolution. We used a probe qubit as 
a system coupled to an environment modeled by a quantum 
map. In this setting there is a straightforward relation between 
some measures of NM behavior and the AFA. The study of the 
time evolution of the NM measure has shown that the stronger 
the chaos of the environment (in the PCM larger A), the more 
Markovian the evolution will be, even if the coupling is strong. 
Furthermore, there are two well defined regimes. For short 
times, there is no dependence with N and the NM is measured 
by revivals in the AFA. In contrast, for large times, the mea- 
sure grows linearly with a slope that vanishes as °= 1 /N. Thus, 
in accordance with 02211 . as N — > °° there can be a remaining 
non vanishing value for non-Markovianity, for a chaotic en- 
vironment. The revivals of the LE were recently related to 
NM behavior ifToll . Here we make a more general approach 
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by allowing the bath to be in a thermal state and expressing 
non-Markovianity in terms of the AFA - a quantity which is 
independent of the set of states over which the average is done. 
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